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(54) MPEG PGR Jitter, frequency ofTset and drift rate measurementa 



(57) A method of measuring MPEG PGR jitter, fre- 
quency offset and drift rate uses a selectable, constant 
measurement bandwidth over non-unllbnn PGR arr1>»l 
times and a variable PGR rate. TTie selectable, constant 
measurement bandwidth Is fixed at a frequency that de- 
marcates betwveen jitter and wander. For each received 
PGR value a PGR intenrai is d^ermlned from the arrival 
time of the cunent PGR value and the anival time of the 



preceding arrival Hme. The PGR values, intervals and 
the fixed bandwidth are input to a set of difference equa- 
tions, derived either from a hybrid digital^aiog PLL 
model where the frequency offset, drift rate and jitter are 
extracted frm vartoua points in the PLL feedbadc con- 
trol system or from a Unear filtering approximation to a 
least mean square (LMS) frequency offset and an LMS 
drift rate estimator off a simple second order lime equar 
tion for the PGR values. 
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Desofpfion 

BACKGROUND OF THE INVENTION 

s [0001] The present invention relates to the quality of compressed video signals, and more particularly to a method 
for measuring Motion Picture Engineering Group (I4PEG) Program Clock Reference (PCR) Jitter, frequency oflset and 
drift rate with a selectable, constant measurement bandwidth over non-unHbrm PCR arrival times and a variable PCR 
rate. 

[0002] Recovering a program system timing dodc (STC) at a decoder side of a digital transmission system, such as 
to a digital video broadcast (DVB) system, is necessary to recreate a video signal. In the Motion Picture Engineering 
Group (MPEG) compression scheme to aOow Isochronous recovery of the STC, explicit Primary Clocic Reference 
(PCR) samples of an encoder's program STC count, such as a 27 MHz dodc count, are sent within a Transport Stream 
CTS). It Is required that the PCR values are both (i) corred at the point of origin and 01) not distorted In the Iransmissbn 
. chain to the point of creating problems in the process of decoding the compressed signals. Measuring the interval 
ts between arrival times of the PCR values, the accuracy of the expected values, the arrival time and resulting Jitter 
accumulated on the PCR values, and the transmitted arrival time In the TS is necessary to assure the confidence of 
decodabllity of the program streams. As Jitter, firequency oflset and frequency drift rate are Important parameters for 
the overall process, a clear definition is needed for what is understood as PCR Jitter and a guidance is needed for its 
measurement method. 

20 [0003] From the specifications set forth In iSO/lEC 13818-1 It Is possible to define a limit masic, as shown in Rg. 3, 
for the frequency deviation from the nominal program STC frequency, 27 Mi-lz. Frequency offset is the differenoe 
between the adual value and the nominal ttequency of the STC at the encoder. For 27 MHz the Dmit is set at +/-^10 
Hz or normalized as +/-30 ppra Frequency rate of change or drift rate is the "speed" at which the frequency of the 
program STC varies with time, i.e.. It isthefirst derivative of the frequency with respect to time or the second derivative 

25 of phase with resped to time. The limit Is set to 75 mHz/sec for the 27 MHz program STC or, nonnaliasd to the STC 
frequency of 27 MHz, as 10 ppm/hr. The phase tderance is set to +/-S00 ns which represents the maximum error of 
a PCR value with resped to the value It should have due to its time position in the TS. 

[0004] The limits for frequency offset and drift rate are Imposed for the STC as it is represented by the values of ttie 
conresponding PCR fields. They indude the effects of the program STC and any possible errors in the PCR calculation. 
30 The limit of 500 ns is not Imposed on the program STC, but on the accuracy representing the PCR values with resped 
to their position In the TS. However the PCR emors are fully equivalent to phase and Jitter errors when the PCR values 
are used at the decoding point to reconstrud the STC. 

[0005] A method for testing program STC itequencyenor, drift rate and PCR Jitter Is described in iSO/lEC 13818-8, 
1996. This test Is carried out for one program in the TS at a time. The procedure uses a plot of PCR ^ues, which 
9$ rspresent the program STC time at which the data was transmitted, i^ainst PCR arrival times, as illustrated In Fig. 1 . 
The an-ival tinries are plotted on the X axis and represent accurate time at the receiver. The PGR values are plotted on 
the Y axis. A linear regression line is fitted to N samples and tiie slope of the line minus one Is the f^uenc^ offM. 
The Jitter is defined as the distance of the PCR values abovert)eiow the line. 

[0006] This method has the problem of deterrnining how large N should be to aocumulate enough samples. If N Is 
40 too targe, low ft-equency timing perturbations, l.e.. wander, cause overiy large measured values of ptter when using ttie 
straight line curve fit rather than a quadratic fit. For ttie Dnear fit ttie slope of this straight line Is used to estimate ttie 
frequency offset. But since the frequency offset often drifts or wanders, the Jitter values may be unduly large due to 
ttiese wander effects. Since ttiese values far exceed ttie Jitter spedfication but are tradced by ttie decoder's phase 
locked loop (PLL). If the drift rate is below limits ttiey are of no concern. Therefbre ttiey should not be in ttie Jitter 
45 measurement and should be separately quantified wtth respedto drift rate. The reverse is also tme of ttie measurement 
of drift rate — small timing variations tiiat are well wittiin ttie Jitter spedfication may, and often do, far exceed ttie drift 
rate spedfication since it is proportional to ttie second derivative of the timing or phase ern:>r. Then^fore ttie drift rate 
measurement should be lowpess filtered to not include timing errors witti pealc variations within acceptable Jitter limits. 
It is preferred to do a quadratic curve fit also to verify ttie compilanoe witti firequency drift rate but. If N Is too small, 
so timing variattons wittiin ttie Jitter limit cause the drift rate esttmate to always exceed spedfication. If tiie PCRs do not 
arrive witti suffident regularity, then ttie regenerated STC n«y jitter or drift. The receiver/decoder may even go out of 
lode 

[0007] The timing efteda visible at a TS Interface pdnt may be modeled as indicated In Rg. 2. This Is an Idealized 
model of an encoder/remultiplexer and communication networtc. The components of ttie model in ttie encoder are 
55 8pedfic to the PCR for a particular program stteam, whereas ttie components In ttie networic relate to ttie entire TS. 
[OOOq The model has a system clock frequency osdilator witti a nominal frequency, F. such as 27 I^Hz. but whose 
actual frequency deviates from this by a function, f^gJfpAY This function depends on ttie time (t) and is spedfic to a 
single program atteam (p). This represents Itequency offeet for tiie nominal frequency which according to the MPEG 
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standard should stay between specified tolemnoes. DrHl rate is the rate of change with time fiff^^iy The drift rate 
also should be below limits set by the MPEG standard. 

[0009] The system dock frequency oscillator drives a PCR counter which generates an idealized PCR count. Np^ 
where (p) refers to the specific program stream and (I) refarB to the bit position in the TS. To this Is added a value from 
s a PCR Inaccuracy source, Mpj, to create the PCR value seen In the TS. Pp^ The simple relationship between these 
values is: 

Pp.i = Np.i*Mpj 

10 

[0010] At the output of the encoder the TS Is being transmitted at a constant bit rate, Rnom- The communicalion 
network beyond the encoder varies the arrival time of MPEG packets from their nominal arrival time by introducing a 
variable delay, JOK whtoh Is the delay in seconds fbr the Ith bit In the TS. The nominal arrival time is calculated from 
R,^ and the bit index in the TS. The equation fbr the actual arrlvd time, A{i). is: 

f5 

A(l) = l/R„oin*D + JCl) 

D is e constant representing the fact that the time origin and bit Index origin are arbitrary and that there is a mean delay 
20 through the communication network. J(l) represents the Jitter In the network delay and its mean value over all time is 
defined to be zero— this sets the value of D. J(l) + Mpj Is the overall PCR jitter which has a maximum error specification 
set In the MPEG standard. 

[001 1] A constant bandwidth method Ibr determining the frequent offset, drffi rate and Jitter Is prefenBd so that all 
measurement s^iems measure the same wlues on the same program PCRs at the stated bandwkJth. Also a constant 
25 bandwidth measurement has more physical meaning in temns of how the recovered STC In the MPEG decoder re- 
sponds to PCR errors since the decoder typically has a constant bandwidth and operates independently of PCR arrival 
Intervals or PCR average rate. 

[0012] What is desired Is a method of measuring MPEG PCR Jitter, frequency of^t and drift rate with a selectable, 
constant measurement bandwidth over non-uniform PCR arrfval times and a variable PCR rate. 

30 

BRIEF SUMMARY OF THE INVENTION 

[0013] Acoordingty the present invention provides a m^hod of measuring MPEG PCR jitter, frequency offset and 
drift rate with a selectable, constant measurement bandwkfth over non-uniform PCR arrivari times and a variable PCR 

3S rate. TTie selected constant bandwidth demarcates jitter from wander and is derived from a limit mask based ixpon the 
standanJs specified for jitter, frequency and drift rate in the MPEG standard. The PCR arrival times are measured with 
a reference clock and are converted to measured PCR Intervate, and the measured PCR intervals, PCR values and 
selected constant bandwidth are input to a set of difference equations to derive the jitter, frequency and drift rate values. 
The set of differenoe equations Is derived etther fmt\ a type 11 phase locked loop or from least mean square (LMS) 

40 linear and quadratic estimators. 

[0014] The objects, advantages and other novel features of the present invention are apparent from the following 
detailed descripUon when read in conJunctk)n with the appended dalms and attached drawing. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

45 

[001 5] Fig. lis a graphical diagram view for estimating frequency oifeet and fitter according to the prior art 
[OOiq Rg. 2 is block diagram view of a modd indicating sources of variaflon in PCR values for a program stream 
In an MPEG Transport Stream. 

[0017] Rg. 3 Isa graphical diagram view of a frequency offset drffll andf tter maskacoording tothe presentinvention. 
30 [001 8] Rg. 4 Is a simplified , general diagram view of a measurament fiHer for measuring frequency offeet, dither and 
jitter according to the present invention. 

[0019] Rg. 5 is a bk>ck diagram view of a phase locked loop (PLL) representaition for measuring frequenqr offeetr 
dither and Jitter according to the present Invention. 

[0020] Rg. 6 Is a bk>ck diagram view of an equivalent circuit to thai of Rg. 5 for deriving a set of difference equations 
S5 according to the present Invention. 

[0021] Rg. 7 is a block diagram view of an equivalent drcult using least mean square first and second order estimators 
ibr deriving a set of difference equations according to the present invention. 

Vmzi Figs. 8a and 8b are comparisons of FIR flnear and quadratic regression fiiters with a dlffierentiator plus lowpass 
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filter end dual differentiators plus lowpasa fittBrs reapecQvely acoonflng to fha praaant invention. 
DETAILED DESCRIPTION OF THE INVENTION 

8 [0023] IWo mettioda are deacribed for measuring frequency offset, frequency drift rate and Jitter of an MPEG System 
Timing Clock (STC) as determined by transmitted PCRs, each method being described as a set of simple difference 
equations suitable for Implementation on a computer TTiey effect constant bandwidth filters operating on norHinlfbrm 
sampled data independent of the PCR rate. Each method, as shown In Fig. 4, has two discreet time inputs: (1) the 
decoded PCR values, PCR(n), for a given program stream f^ the MPEG T^spoit Stream (TS); and (2) the PCR 

to inten^als, PCR.int(n), between PCR values for the given program stream aa measured wHh an asynchronous precision 
refisrence oscillator from the PCR arrival times. Each method has a single control parameter — the measurement 
bandwidth (BW) which is a demarcation frequency between Jitter and wander - that Is selectable and constant inde- 
pendently of the average PCR rate. Each n>ethod has three discrete time output measurement results: (1) STC fre- 
quency ofliiet, FO(n); (2) STC finequency drift rate, DR(n): and (3) PCR Jitter, Jit(n). 

10 [0024] The equations for the first method are derived from a hybrid digital-analog phase lodged loop (PLL) model, 
shown in Rg. 5, where the three output paremetere are extracted from various points in the resulting faedbad( control 
system. Additional high-pass tittering and scaling Is added to Jlt(n) to Improve separation of wander components and 
scale the output to units of time per second. Additional low-pass filtering Is added to FO(n) and DR(n) to Improve 
rejection of jitter components. A practical set of difference equationa are derived from the PLL model such that frequency 

20 of^et, drfft rate and Jitter may be measured with a setodabie, constant bandwidth independently from the non-uniform 
PCR an-lval times. 

[0025] The type II PLL shown in Rg. 5 has two ideal Integrators - a digital accumulator and a voltage controlled 
oscillator. TTiis creates a second order highpass dosed loop response at the output of the phase subtraction. Therefore 
below the loop bandwidth the response is proportional to drill rate and at»ove the loop bandwidth the response Is 

25 proportional to jitter. A first order highpass filter is added to the Jitter measurement to remove the efiects of drift rate. 
Conversely a lowpasa filter Is added at the drift rate output to remove the efTects of Jitter ftom that measurement. TTie 
filters are best Implemented using digital signal processing (DSP) on the raw data and. since the PCR rate Is the sample 
rate, the average PCR rate is detennined and the confidents adjusted to maintain approximately the correct meas- 
urement filter bandwidths and control loop damping. It Is possible to convert the PLL to a fiiily sampled system drWen 

30 by a measured PCR interval using a firee running oscillator and the decoded PCR values such that the bandwidth is 
Independent of non-uniform PCR rate and the measurementvaluesarecomputed In software without an expiidtvoltage 
controlled oscillator. 

[0026] The VCXO may be constructed from a subsystem having a free^ning OCXO of 0.1 ppm accuracy and a 
FLL lodclng a VCXO with +^30 ppm renga This may be used to calibrate the frequency offiset output to 0.1 ppm 
35 accuracy If desired. Otherwise the VCXO may be used alone and If s frequency error or offset verified by applying a 
known, accurate frequency, and subtracting the error ftom subsequent measurements. This desaibes a PLL imple- 
mentation as a hybrid DSP and analog signal processing method. 

[0027] A practical set of diflerenoe equations may be derived from the PLL model duplicating the closed-loop re- 
sponse such that frequency offset, drift rate and jitter may be measured with a selectable, fixed bandwidth Independently 

40 of the non-uniform PCR arrtval times. Such an equivalent circuit is shown In Rg. 6 for deriving such equations. TTie 
PCR(n) values are input to a differentiator, and the output from the diflerentlator Is Input to a second order lowpass 
filter for non-unifonn sample intervals. Also Input to the second order lowpass filter is the cone^ndlng PCR.Int(n) 
from the last received PCR(n-1 ). The filter output y(n) Is multiplied by the overall open loop d.c. gain. Godc, to obtain 
a normalized frequency from which unity is sutstraded to derive the bi-polar frequency offset FO. The filter outputs for 

45 two consecutive samples, y(n) and y(n-1), are Input to a second differentiator to produce a jitter error output J(n), 
which is in tum input to both a first order lowpass filter to derive the drift rate. DR(n). and a first order highpass filter to 
derive the jitter, Jlt(n). From this iiiustrated drouit the dlfierence equations may be derived: 

^ y{n)=2*y(n-1)-y{n-2)-(y(n-1)-N1Mn-2))*6.932*BVV*PCR«^ 

where Fo Is 27MHz, BW Is the measurement bandwidth in Hertz and N1 Is a constant equal to 1/(1 +pi*BW/60). From 
y(n) frequency offset. FO(n), is y(n)*Godo-1 where Qodc Is a constant equal to (1-Nir5.932*BW. Jitter, Jit(n), becomes 
Jlt(n}=A1*(Jlt(n-l)+J(nH(n-1)) where J(n)=y(n)-y(n-1) and A1 Is exp(-2*pl*BVTO.68/30). Rnaliy drift rate, DR(n), Is 
65 AT DR(n-1)4-(Qodcn'prJ(n) where Tp Is 1/(2>rQ1) and A2 Is exp(-1/30*Tp). In this way DR(n) Is post-filtered with a 
fixed bandwidth of 0.01 Hz and y(n) Is Ntialteed to 1/Godc to minimize settling time. 

[0028] The second method Is derived from a linear filtering approximation to a least mean square (LMS) frequency 
offset estimator (Ilnear-temn coefficient) and LMS drift rate estimator (quadratic4erm coefficient) of a simple second 
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order HmQ equation for the PCR values. The resulting difierenoe equations yield results simllar to those in the method 
described above, but are a little simpler to compute. 

[0029] Given PCR.AT(n) as the measured arrival time of each of the received PCR(n) values using an accurate 
reference clocic, the following eqi^tJon defines the expected PCR(n) values with terms for the program decoder STC 
0 fhsquency offsett drift rate and Jltten 

PCR(n)==a♦PCf^T(n)•F8tc+<PCF^T(n))^2^^^e(PCR^AT(n)) 

10 cycles of Fsta The last term Is the Jitter as a function of time. PCR_AT(n) Is a random variable added to PCR(n) due 
to the comt>ined effects of anival time Jitter and PCR(n) encoding error. It to assumed to be a zero-mean, Independent, 
Identically Distributed random variable, which Is falriy realistic. Therefore the PCR(n) values are themselves I.I.D, 
random variables, a Is an Initial phase error relative to the reference osdilator, Fate is the system timing clock frequency 
and c Is one-half the drift rate since drift rate (IHz/s) is a second time derivative of phase, DRB2?e. The known ideal 

15 system docic frequency. Fo^MHz, is subtracted to find the STC frequency offset and the terms are normalized by 
dividing by Fd. Since the STC drift rate is typically very small relative to the allowable frequency offset values and Is 
preferred to be included within the frequency offset measurement, the equation may be simplified for analysis of fre- 
quency offset, b, without the drift rate term, c. The drift rate and higher order derivatives affect the frequency offeet 
value rather than defining frequency offset with higher order components removed. 

20 

PCK_erT(n)=a+PCPOKT(n)*b+€(PCRJkT(n)) seconds 

where b=(Fsto-Fo}/Fo. The Jitter term now is in seconds of Jitter and a is the initial phase offeet in seconds and b is the 

25 frequency offset of the STC as a dimenstoniess ratio. This equation, negteding drift rate, Is the equation for F^R(n}8 
allowing a linear regression or straight line fit of the data to be determined from estimates of the coefficients a and b 
to minimize the mean square error due to Jitter. This provides the dual result of first detemilning STC frequency offeet 
FO as the LMS estimate for b, and second the Jitter by subtrecting the FO effect due to b and the Initial phase value 
a as suggested by the IMPEG standard. 

30 [0030] in order to determine the number of PCR(n) values with which to compute the straight line fit, a comparison 
of the linear regression analysis to a linear filtering process fe necessary. This way the effective bandwidth of the 
averaging of the linear fit to the PCR(n) data may be found fer an averege PGR rate where the PCR(n) values are 
temporarily assumed to arrive at a uniform average rate. This uniform rate assumption is Justified In this context due 
to the large number of samples that are used in order to properly estimate FO within the oonstrainds imposed by the 

^ very small drift rate tolerance of the MPEG standard, 10 ppm/Hr. Once a good approximation of a filter response to 
the straight line fit regression analysis is determined, a constant t>andwidth filtering process may be constructed for 
non-uniform sampled data that is itself a function of the instantaneous sampfe interval In order to develop a constant 
bandwidth measurement that is Independent of PCR rate. This is important since it allows a continuously updated 
filtering process to be mechanized to find practical and consistent estimates of the PCR Jitter, STC fluency offset 

40 and drift rate for comparison to the published limits of the MPEG standard. Further since the bandwidth of the meas- 
urement is constant, a standardized demarcation frequency based on the STC drift rate and PCR accuracy specifica- 
tions may be used to determine a practical and standardized filter t>andwidth that has a direct relationship to the re- 
covered STC clodc characteristics In the MPEQ decoder. 

46 UnearReormsion AnaiyBia: 

[0031] Ijet TO be the average PCR.Interval determined by the reference dock and assume uniform arrival times 
spaced TO seconds apart. 

so 

PCK.err(n)=a+b*TO*n+e (n*TO) seconds 
Given N samples of PCR values, least mean square (LMS) estimates of a and b are given by vector B: 

B=[a;b]; B=r^*(N*NV*(N^*X) 

where 
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T"''=[1,0;OJO''*1; (hrNY=I2(2N+1),-6;^,12/(N*1)ni/(N(N-1))) 

(N^*X):= [2^,^PCR.err(n); D^^^PCR^efrCnrn]. 
Initial phase la determined below: 

a=2/(N*(N-1)T0){r„„i4,PCR.err(n)*(2N*l^)}. 
Frequent^ offeet estimate for N uniformly spaced samples: 

bs 6/(N*(N-1)TO){X:^^i^PCFLen<nn2n/(N+1)-1)}. 
Quadratic Regrasslon Analysis: 

[0032] Without neglecting the STC frequency drift rate find an 1^ quadratic curve fit and estimate of the coefflclent 
c. where c=(STC frequency drift rBte)/2sDR/2: 

PCR.err(n)=a+b*TO*n+c*rO^*n^+e{n*TO) seconds 

B=[a; b; c]=r^N*N V(N^*X); r^1.0,0; 0 J0'\o; OATO"^; 

(M*n'^)"^=P(3(N+1)+2),-18(2N+1)^ 

-18(2N+1),{12(2N+1)(8N+11)l/{(N+1XN+2)K-180/(N+2); 

30,-180/(N+2), 180/{(N+1)(N+2)}]/{N(r4-1Hr^2)}; 

(N'^*X)=[X^i^PCff^err<n); 2^^i^PCR^erT(nrn; 2:^^PCl^err(nrn^. 

TTien LMS estimate for frequency drift rate for N uniformly spaced samples is: 

DR=60/{N(N-1)(N-2}TO^}*2n^^PCR^err(n){6n^/{(i^1)(N+2)}-6n/(N*2)+l}. 

[0039] These two LMS estimators, frequency offset from linear regression and frequency drift rate from quadratic 
regression, may be implemented as finite impulse response (FIR) filters. First describing the LIMS estimate of the 
frequency offset In terms of a FIR filter as a function of z and itie number of PGR values averaged, N: 

Hr<z.N)=6/{N(N-1)TP}*2"^*Lfc^fi.iyZ*'{2(lc+l)/(N*)-1}. 

TTie number of samples used to estimate FO Is iM=Pr/(2*BW} where Pr is an average PGR rate, TO is 1/Fs and BW is 
the measurement bandwidth. Therefore as an example N=floor(Pr/(2*BVV)) = 125 samples with Pi=24 and BW=0.1 Hz. 
Second describing a FIR filter as a function of z the frequency drift rate estimate Is: 

Hdr(z.N)K6(V(N(N-1)(ri-2)TO^)}z"^Zk^i{(6(lc*1)^«^ 
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[0034] For comparison a second order recursive Butterworfh lowpass filter of bandwidth BW is cascaded with a 
differentiator. A cascade of the recursive LPF and the differentiator response Is compared to the FIR LMS frequency 
offset estimator in Fig. 8a. Drift rate RR response Is shown In Rg. 8b where: 

N=floor(Pr/(2''°*Bw)=nuirtber of samples 
processed by the FIR filter at the average PGR rate of Pr. 

For comparison a 3"^ order recursive Buttenvorth LPF of bandwidth BW cascaded with two difTerentiators is compared 

10 to the FIR LMS drift rate estimator In Fig. 8b. 

[0035] From this analysis there Is a close approidmatlon of linear and quadratic regression of uniformly sampled data 
to two linear, recursive littering processes cascaded with simple differenliatDrs. Given an estimate of the average PGR 
rate, Fr, the bandwidth, BW, of the equivalent recursive filter process to replace a Unear regression analysis of N 
samples Is BWsPr/(2*N) and lor quadratic regression BW=Pr/(2'"2*N). These equivalent processes may be constructed 

18 from difference equations for differentiators and ButtenArorth LPF filters as shown above. Since the equivalent recursive 
filters only require storage of two prevtous PGR values and relatively minor multiply/adds, they are much easier to 
compute and may be done concurrently for many program PGRs with a single prcx^ssoryieidlng a constant t>andwldth 
measurement result The PGR jitter may be detemnined by subtracting the estimated PGR values due to fluency 
offset by using the estimated value from the filtering process and applying a HPF to reject the drift rate components 

20 from the result. 

[0035] An equivalent circuit for this second method la shown in Rg. 7. The PGR^Int(n) is derived ftom the PGR^T 
(n) times and subtracted from the differentiation of the PGR(n) values, nomnalized, to produce a PCR(n) en^or value, 
Foo(n). Foo(n) Is Input to a first Onear LMS estimator which includes a derivative function and a second ord^ lowpass 
filter to produce the frequency offset, FO(n). The frequency offset Is Input to a second quadratic LMS estimator which 
25 includes a derivative function, differentiator and first order lowpass filter to produce the drift rate. DR(n). The frequency 
offeet also Is multiplied by the PCR^Int(n) and the result subtracted from the PCR(n) enor value to produce an instan- 
taneous Jitter value, JJ(n). The Instantaneous Jitter value is then processed by a cascade filter, having first a first order 
lowpass filter fbllowed by afirst order Mghpass filter, that produces the PGR jitter, Jit(n). The resulting difference equa- 
tions are: 

30 

FO(n)=PCR^Int(nnn*BW)^Too(n)*E(n)*2E(nrFO(n-1VE(n)^T^^ 

where 

35 

Foo{n)=IPCR(n)-PCR(n-1B/Fo-PCR^int(n); 

40 E(n)=expI-PGR^int(nrn-'BWl (exp(-x)-1 W/2-x^/3): 

DR(n)=E(nrDR(n-1>+n*BWTF0(n)-F0(n-1)l; 



JJ(n)=Foo(n)-FO(n)*PCR^Int(n): 



Jit(n)=E(n)*[JJ(n>JJ(n-1)*2*Jil(n-1)-E(nrJlt(n-2)]. 

so 

[0037] Thus the present invention provides a method of detemilning MPEG PGR Jitter, fluency offset and drift rate 
ftom the PGR values arriving at non-unltonn Intervals and a selecteble, constant bandwidth which Is the frequency that 
separates wander from Jitter using a set of difference equations derived either ftom a type II phase locked loop or from 
LMS first and second order estimators. 
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Claims 

1 A method of measuring program dock reference (PGR) Jitter, frequency offset and drift rate with a selectable 
constant measurement bandwidth over non-unSfbrni PGR arrival times and a varfal>l8 PGR rata comprising the 
steps of: 

measuring an interval between PGR arrfvat times with an asynchronous precision reference osclDaton 
fixing the selectable constarit measurenoent bandwidth Independently of the variable PGR rate, the fixed se- 
lectable constant measurement bandwidth demarcating a frequency between Jitter and wander; and 
calculating the PGR jitter, frequency offeet and drift rate as a function of the i nterval between PGR arrival times, 
the PGR values at each PGR arrival time and the fixed selectable constant measurement bandwidth In the 
form of a set of difference equations. 

Z The method as recited in claim 1 wherein the set of dlfferenoe equations Is derived from a hybrid digital-analog 
phase lodced loop (Pi±) modei where the jitter, frequency offset and drift rate are extracted from various points in 
the feedbadc control system for the PLL 

3. The method as recited in claim 1 wherein the set of difference equations is derived from a Onear filtering approx- 
imation to a least mean square (LMS) frequency ofliBet and an LMS drift rate estimator of a simple second order 
time equation for the PGR values. 
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